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Abstract 

Presently the quarkonium spectrum, written in terms of the quark MS-mass, is known 
at 0{o? s m) accuracy. We point out that in order to achieve O(o.gm) accuracy it is 
sufficient to calculate further (I) the 0(a s m) relation between the MS-mass and the 
pole-mass, and (II) the binding energy at 0(a s m) in the large-/?o approximation. We 
calculate the latter correction analytically for the lS-state and study its phenomenological 
implications. 
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1 Introduction 



In recent years theory necessary for precise description of a heavy quarkonium such as bot- 
tomonium or (remnant of) toponium has developed significantly. In particular, development in 
the computational technology of higher order corrections to the quarkonium energy spectrum 
[0], and subsequent discovery of the renormalon cancellation || [| enabled accurate determi- 
nations of the MS-mass of the bottom quark jjj |5], 0, || and (in the future) of the top quark 
@. In these determinations a major part is played by the spectrum (mass) of the lS-state 
quarkonium.f] 

It is legitimate to consider that the present perturbative calculation of the quarkonium 
spectrum, when expressed in terms of the quark MS-mass, has a genuine accuracy at 0{a% m). 
In fact, in formal power countings, the last known term in the relation between the MS- 
mass and the pole- mass of a quark is (9(a|m), while the last known term of the binding energy 
(measured from twice of the quark pole-mass) is 0{a i s m). The former term includes in addition 
to a genuine 0{a\m) part the leading renormalon contribution which does not become smaller 
than O(Aqcd) [0- This renormalon contribution is cancelled ]3], |J against the renormalon 
contribution |T1J contained in the latter term. Therefore, after cancellation of the leading 
renormalons, the genuine 0{a 3 s m) part of the mass relation determines the accuracy of the 
present perturbation series relating the quark MS-mass and the quarkonium spectrum. 

As stated, for the binding energy the calculation including the genuine 0(a A s m) correc- 



tions has already been completed. Also, the "large-/?o approximation" [|12| is known to be a 
pragmatically feasible and empirically successful estimation method of the leading renormalon 
contributions. Taking these into account, one finds that it is sufficient^] to calculate further the 
following two corrections in order to improve the accuracy of the spectrum by one order and 
to achieve a genuine accuracy at 0{a^m): (I) the (9(a|m) relation between the MS-mass and 
the pole-mass, and (II) the binding energy at (9(a|m) in the large-/?o approximation. This 
is because the leading renormalon contribution in the full 0(a 5 s m) correction to the binding 
energy will be incorporated by the large-/9o approximation and the remaining part is expected 
to be irrelevant at 0{a A s m). 

Of these two corrections we calculate (II) analytically for the lS'-state in this paper. Then 
we study its phenomenological applications. We can check validity of the above general argu- 
ment explicitly at (9(a|m) where we know the exact result. This will also be demonstrated. 

Written in terms of the quark pole- mass m po i e , the mass of the quarkonium lS'-state is 
given as a series expansion in the MS coupling constant «s(/i) defined in the theory with ni 
massless quarks: 

M 15 = 2m pole - l as (//) 2 m pole ^(^^) n p n (L), (1) 

n=0 



* See e.g. || || for introductory reviews of the subject. 

t There exist other methods in which the renormalon contribution contained in the pole-mass is subtracted in 
certain approximations (e.g. [0). In our opinion our method is a most natural one, embedded in the algorithm 
for higher order calculations of the whole spectrum; cancellation of infrared sensitivities in the whole spectrum 
follows from the fact that the quarkonium is a color-singlet small-size system. 
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where P n (L) is an n-th-degree polynomial of L = log[3/x/(4as(/i)m po ic)]- At each order of the 
perturbative expansion c n = P n (0) represents a non-trivial correction, while the coefficients of 
L's are determined by the renormalization-group equation.^ The polynomials relevant to our 
analysis read 

(2) 
(3) 

+ \Pi + l^ocA L + c 2 , (4) 

+ Q/?o - \PoPi + ^2 - \ftc x + ^icx + 2f3 c^j L + c 3 . (5) 

/3 n 's denote the coefficients of the QCD beta function given by 

2 n 38 2857 5033 325 2 
fa = \\--n h ^ = 102--^, fo = — 2 W ni + ^A 1 ' () 

Note that, in the renormalization-group evolution, running of the coupling included in log as 
should be taken into account properly. This feature is unique to the perturbation series of a 
nonrelativistic boundstate spectrum. Presently Cfi S £1X6 known up to n = 2 0, : 

97 11 

1793 2917tt 2 9vr 4 275C 3 / 1693 llvr 2 19C 3 \ / 77 tt 2 2( 3 \ 2 

Co = 1 1 h )ri] + 1 1 nf. 

12 216 32 4 V 72 18 2 / V108 54 9 / 1 

(8) 

The aim of this paper is to calculate c 3 in the large-/?o approximation. 

One might think that alternatively c 3 may be estimated using the asymptotic form of 
the series expansion of the QCD potential at large orders. We could not, however, find a 
justification that c 3 can be estimated with 0(a s m) accuracy in this manner. 



2 Outline of the calculation 

In the large-/?o approximation the wave functions and energy spectra (measured from 2m po i e ) 
of quarkonium states are determined by solving the nonrelativistic Schrodinger equation 

ip{x) = Etp(x). (9) 

Here, m denotes the mass of the quark and it is irrelevant whether we use the pole-mass 
or the MS-mass for m within our approximation; in particular it does not affect the leading 

For n > 3, c„'s include powers of log as unrelated to the renormalization group, i.e. which is not accom- 
panied by log . 
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Table 1: First few coefficients of the perturbative expansion of R(r). Here, t = log(/xr) + | + 7_b- 



renormalon cancellation. Vg denotes the QCD potential in the large-/?o approximation given 



as a perturbation series [11 



V Po (r) = -C F ^-xR(r) 



R(r) = £ 



n=0 



/ /3 ag(/i) \ w 

v 4vr ; 



(10) 



CV = 4/3 is the color factor. First few coefficients R n (fir) are listed in Table [I]. To obtain 
S- wave solutions, we introduce a dimensionless variable z = Cpctsmr and set 



ip = - exp 

z 



-J dzW(z) 



E = {CpasYme. 



Then the equation becomes 

W - W 2 - R/z = e. 
We expand W, R and e in perturbation series as in ( |i~0"D and 

'i3 a s 



n=0 



n=0 

and substitute them to (|i~2"l) . We find 

n 

For the lS'-state the zeroth-order solution is given by 



47T 



w 



1 1 



2 2' 

and we may solve ( |T4] ) recursively for n > 1: 



Eq = — 



WJz) 



-I l>00 I 

-- \ dU 2 e- t [Y^W k {t)W n ^ k {t) + Rnft/t 

z poo /n~* 
-- / dttV Q]^(t)^_ fc (t) + iJ n (t)/H£ 



(12) 



(13) 



(14) 



(15) 



(16) 
(17) 



Thus, e n 's can be obtained by evaluating multiple integrals. We may render these integrals 
to forms which resemble Feynman-parameter integrals that appear in calculations of multi- 
loop Feynman diagrams. In this way we could use various techniques developed for Feynman 
diagram calculations.^ We obtained the first two terms of the perturbative expansion as 



e 2 



-(L + l) 



-(3L 2 + 4L + 1 + — + 20 
6 



(18) 
(19) 



where L = L + | = \og[fi / [C FOtsm)} + |. From these we confirmed the corresponding parts of 
Ci and C2 in eqs. ([7] 

^3 = - 

It follows that 

c 3 (large-/? ) = 0$ 



. We also obtained a new result: 
LU- I:'/. ( — + I6C3J L - 2 + tt 2 + I6C3 + — - 2vr 2 C 3 + 24( 5 



, , 517 , 19 vr 2 , 11 C 3 T 



864 + 144 
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4 k 2 Cs | 3C 5 



1440 



(20) 



(21) 



We note two points: (i) The result includes the level-5 zeta values (£5 and tt 2 ^). (ii) The result 
does not include je, that is, all 7^ which appeared at intermediate stages of the calculation 
got cancelled. A more detailed description of our calculation will be published elsewhere. 



3 Phenomenological applications 

We examine the series expansion of the quarkonium lS'-state spectrum numerically for the 
bottomonium and (remnant of) toponium. The input value for the coupling defined in the 
theory with 5 massless flavors (ni = 5) is ap(Mz) = 0.119. We evolve the coupling and match 
it to the coupling of the theory with ni = 4 following \TE\. Also we take the input values of 
the MS-mass, m = ^Ms( m Ms)' an< ^ the pole-mass as fn b = 4.20 GeV/m bjPO i c = 4.97 GeV and 
m t = 165 GeV/m 4jPO ie = 174.79 GeV. The natural size of a quarkonium system is the Bohr 
radius. We define a corresponding scale parameter \i B such that 

fx B = C F as(n B )m pole (22) 

holds. We examine the series expansion ([!]) with two different choices of the scale \x in Table |2|. 
The last terms (with stars) are evaluated using the value of C3 in the large-/?o approximation 



pi). One sees that for the bottomonium the series expansions do not converge at all. For the 
toponium the series expansions converge very slowly. The numbers in square brackets represent 
estimates of the last terms using the asymptotic form of the series expansion of the potential 



- C F ^ x R n ^r) ~ - 2 e5/6g ^M H x 2^n\ for n » 1. (23) 
r 7r 



We find the techniques developed in |tj particularly useful. 
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Series expansion of Mi 5 [GeV] 





A* = Mb 

[ expansion parameter = as(^s) ] 


H = m 

[ expansion parameter = as(rn) ] 


Bottomonium 

(ni = 4) 


2 x (4.97 - 0.11 - 0.12 - 0.20 - 0.25*) 

[-0.26] 


2 x (4.97 - 0.06 - 0.08 - 0.11 - 0.15*) 

[-0.16] 


Toponium 

(ni = 5) 


2 x (174.79 - 0.77 - 0.35 - 0.25 - 0.13*) 

[-0.14] 


2 x (174.79 - 0.46 - 0.39 - 0.28 - 0.19*) 

[-0.25] 



Table 2: Numerical evaluation of the series expansion of Mis, e Q- ([!)• Relevant parameters for 
the bottomonium are: m& = 4.20 GeV, m& )PO i e = 4.97 GeV, hb = 2.05 GeV, Og (^b) = 0.309, 
a s(™b) — 0.229; those for the toponium are: m t = 165 GeV, mt iPO i c = 174.79 GeV, \i& = 32.9 GeV, 
a { g\n B ) = 0.1411, af{m t ) = 0.1092. 



The series approaches its asymptotic form faster when we choose \x = \xb- These features are 
consistent with dominance of the leading renormalon contributions. 

Next we rewrite the series expansion of Mis m terms of the MS-mass instead of the pole- 
mass. The leading renormalon contributions cancel in this case. Presently the relation between 
the MS-mass and the pole-mass is known up to three loops [fL6[ , |17|] . The only scale in this 
relation is the quark mass. Thus, we have two choices of scales, m and /is, in writing the series 
expansion of Mis; we take // = m below. We eliminate the pole- mass completely and expand 
in as(jn). We should properly take into account the fact that the renormalon contributions 
cancel between the terms whose orders in as differ by one [18]. To this end we proceed as 
follows. We rewrite (|1|) as 



Mis = 2m po i c x F 1 + ^ Qn a,s{m) n 



n=2 



2m x 



1 + J2 dn 



a s (m 



n=l 



1 + Qnas(m) n 



n=2 



(24) 



where Q n 's are polynomials of log [as (m)] and <i n 's (the n-loop coefficients of the mass relation) 
are just constants independent of as(m). We identified Q n a s as order a 
the last line to a single series in as- Numerically we findT] 



g 1 and then reduced 



Mis = 2 x (4.20 + 0.36 + 0.13 + 0.040 + 0.0051 fl ) GeV (Bottomonium), (25) 
Mis = 2 x (165.00 + 7.21 + 1.24 + 0.22 + 0.052 fl ) GeV (Toponium). (26) 



The last terms (with sharps) are evaluated using the values of C3 and [12| in the large- 
Po approximation. Convergences of the series improve markedly in comparison to those in 
Table |2|. (Previous analyses similar to the one presented above can be found in [f| 0] and 
references therein.) 

As we argued in Section [l], parametric accuracy of the last terms in (251,261) is (D(a s m) 
and we need to know further only the exact value of to make a perturbative evaluation 

^ The values of d\ ~ d 3 are taken from eq. (14) of [p7| . 
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accurate up to this order (the exact form of C3 is not necessary). In order to verify validity 
of this argument, we replace C2 by its value in the large-/3o approximation. Then the 0{oP s m) 
terms of fl2"5] , |2"6l ) change to 0.043 and 0.22, respectively. Thus, we do not lose accuracy at this 
order by the replacement. On the other hand, if we replace C2 and ofa by their values in the 
large-/9o approximation, the same terms change to 0.056 and 0.31, respectively. Thus, we lose 
the accuracy at C(a|m). These aspects are consistent with our general argument. Also, they 



suggest that the last terms of (|25| , |26|) would be reasonable estimates of the orders of magnitude 
of the exact O^a^m) terms. 

Finally we examine if we can use the asymptotic form of Vp (r), eq. (p3|), to estimate 
C3 (large- /?o). From the asymptotic values of the last terms for /z = in Table |^, we may 
extract approximate values of c 3 (large-/3 ) as 



■ 2.55xl0 3 (n, = 4) 
c 3 (asympt)H ^ gg x xq 3 (n, = 5) ' (27) 



while the corresponding values of C3(large-/?o) are, respectively, 2.46 x 10 3 and 1.91 x 10 3 . If 
we substitute C3(asympt), the last terms of ( |25|J26| ) change to 0.0034 and 0.051, respectively. 
Therefore, we see that in the case of bottomonium the use of the asymptotic form does not 



reproduce our result in the large-/3o approximation (|25f ) with good (relative) accuracy. Presum- 
ably it is a sign of the next-to-leading renormalon contribution, which does not become smaller 
than ~ Aq CD ■ (Aqcd//^) 2 ~ 2 MeV, and which is not included in c 3 (asympt). Namely, we 
conjecture that already the last term in eq. ( ^5|) stands close to the limit where an improvement 
of convergence of the perturbation series is possible by cancellation of the leading renormalon 
contributions. 



4 Conclusions and discussions 

We have calculated the 0(a 5 s m) correction to the quarkonium lS-state energy spectrum ana- 
lytically in the large-/3o approximation. As a result, in order to predict the 15-state spectrum 
at 0{a\m) perturbative accuracy, only the four- loop relation between the MS-mass and the 
pole-mass remains to be computed. Within the present approximation, the perturbation series 
of the bottomonium and toponium spectra show healthy convergent behaviors up to a genuine 
0{a%m) accuracy. 

In the case of bottomonium, current theoretical uncertainty due to non-perturbative effects 
is estimated to be of the order of 0.1 GeV || §. It is much larger than the size of the last 
term of (|25|). We hope that in the future the theoretical uncertainty due to non-perturbative 
effects will be reduced by applications of e.g. lattice calculations or combinations of operator- 
product-expansion and sum rules. 

Part of the genuine 0{a\m) corrections have already been calculated [EJ. Taking // = m, 
their individual sizes are evaluated to be AM15 ~ ±0.05 GeV both for the bottomonium and 
toponium; if we take \i = they become even an order of magnitude larger. We do not 
know as yet whether it may indicate a breakdown of perturbative expansion of the spectrum 
at 0(a 5 s m) or a requisiteness for some new cancellation mechanism. It is also possible that 
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the sum of all the genuine 0(a 5 s m) corrections turns out to be much smaller. In any case, 
it would be better to separate the discussion of this problem from the determination of the 
genuine 0{a i s m) corrections (as we advocated in this paper), rather than to regard them as 
inseparable constituents of C3. 

Acknowledgements 

Y.K. was supported by the Japan Society for the Promotion of Science. Y.S. was supported 
in part by the Japan-German Cooperative Science Promotion Program. 

References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7 

[8] 
[9] 

[10 

[11 
[12 

[13 

[14 
[15 



A. Pineda and F. Yndurain, Phys. Rev. D58, 094022 (1998). 

K. Melnikov and A. Yelkhovsky, Phys. Rev. D59, 114009 (1999). 

A. Hoang, M. Smith, T. Stelzer and S. Willenbrock, Phys. Rev. D59, 114014 (1999). 

M. Beneke, Phys. Lett. B434, 115 (1998). 

A. Penin and A. Pivovarov, Phys. Lett. B435, 413 (1998); Nucl. Phys. B549, 217 (1999). 

A. Hoang, Phys. Rev. D61, 034005 (2000); M. Beneke and A. Signer, Phys. Lett. B471 
233, (1999). 

A. Hoang, et al., "Top-Antitop Pair Production Close to Threshold: Synopsis of Recent 
NNLO Results", |hep-ph/0001286| . 



M. Beneke, "Perturbative heavy quark- antiquark systems", [hep-ph/ 99 11490 



Y. Sumino, "Renormalon Cancellation in Heavy Quarkonia and Determination of m^, 
m t ", |hep-ph/0004087. 



M. Beneke and V. Braun, Nucl. Phys. B426, 301 (1994); I. Bigi, M. Shifman, N. Uraltsev 
and A. Vainshtein, Phys. Rev. D50, 2234 (1994). 

U. Aglietti and Z. Ligeti, Phys. Lett. B364, 75 (1995). 

M. Beneke and V. Braun, Phys. Lett. B348, 513 (1995). 

N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Lett. B470, 215 (1999); B. Kniehl 
and A. Penin, Nucl. Phys. B563, 200 (1999); Nucl. Phys. B577, 197 (2000). 



D. Broadhurst, frep-th/ 9604123$ Eur. Phys. J. C8, 311 (1999). 



K. Chetyrkin, B. Kniehl and M. Steinhauser, Phys. Rev. Lett. 79, 2184 (1997). 



8 



[16] K. Chetyrkin and M. Steinhauser, Phys. Rev. Lett. 83, 4001 (1999); Nucl. Phys. B573, 
617 (2000). 

[17] K. Melnikov and T. v. Ritbergen, Phys. Lett. B482, 99 (2000). 

[18] A. Hoang, Z. Ligeti and A. Manohar, Phys. Rev. Lett. 82, 277 (1999); Phys. Rev. D59, 
074017 (1999). 



9 



